Effects of metallic contacts on electron transport through graphene 
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We report on a first-principles study of the conductance through graphene suspended between Al 
contacts as a function of junction length, width, and orientation. The charge transfer at the leads 
and into the freestanding section gives rise to an electron-hole asymmetry in the conductance and 
in sufficiently long junctions induces two conductance minima at the energies of the Dirac points for 
suspended and clamped regions, respectively. We obtain the potential profile along a junction caused 
by doping and provide parameters for effective model calculations of the junction conductance with 
weakly interacting metallic leads. 
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Introduction. Graphene is a two-dimensional allotrope of 
carbon with atoms arranged in a honeycomb lattice. Suc- 
cessful fabrication of graphene monolayers by means of 
mechanical exfoliation of graphite [ij or epitaxial growth 
on silicon carbide Q has ignited tremendous interest in 
this material Its favorable characteristics such as 
chemical inertness, low dimensionality, extremely high 
mobility, and easy control of carriers by applied gate 
voltages, along with patterning using nanolithography, 
open possibilities for further miniaturization of devices 
and the emergence of a carbon-based "post-silicon" elec- 
tronics 0,11]. 

The electron transport at the nanometer scale is signifi- 
cantly affected by the contacts. The role of metallic leads 
in determining the transport properties of graphene- 
basedj unctions has been addressed by several theoret- 
ical mm and experimental [lll-[l5j| studies. Yet, a 
parameter-free description of transport in these systems 
is still lacking. Previously, the leads were modeled as 
infinitely doped graphene regions that support a large 
number of propagating modes In this model, the im- 
portant parameter is the ratio W/L between the width 
(W) and length (L) of a junction: at the Dirac point, the 
universal, diffusive- like regime is reached for W/L ^ 1. 
Additional studies within a tight-binding approach have 
been reported, in which either the leads form a square 
lattice [a, 0] , or the coupling to the leads is modeled via 
energy-level broadening Q. 

In this Letter, we report on large-scale first-principles 
calculations of electron transport in suspended graphene 
taking into account the effects of metallic leads. We study 
a non-magnetic junction made of graphene contacted un- 
derneath by two aluminum (Al) leads with a small in- 
plane mismatch of less than 1%. This junction is a pro- 
totype system with weak interaction between the metal 
contact and graphene, with no covalent bonds formed. 
Our work represents the first quantitative study of elec- 
tron transport through metal-graphene junctions to ex- 
amine in detail previous models [1, [sj with dimensions 
relevant to experiment lll-ll4|. up to 100 nm wide and 



13.6 nm long with various W/L ratios. 

The difference in the work functions of the metal and 
graphene leads to the charge transfer and doping of the 
graphene layer [16|. We obtain the corresponding po- 



tential profile generated by the doping of the graphene 
layer along the junction. Not surprisingly, the potential 
in the region contacted by the metal starts deviating from 
its bulk value before reaching the geometrical edge. Fi- 
nite doping results in specific transport features for long 
junctions: two conductance minima appear at the ener- 
gies of the Dirac points of graphene in the clamped and 
suspended regions, respectively. For shorter junctions, 
where the two minima cannot be resolved, an electron- 
hole asymmetry in the conductance is still appreciable. 
The impact on transport of the orientation of graphene 
in the junctions is found to be negligible away from the 
Dirac point. Using the potential profile obtained from 
our first-principles calculations and a small energy broad- 
ening in the self-energies at the leads, we demonstrate 
that a TT-electron tight-binding model can accurately re- 
produce our first-principles transport results. This en- 
ables us to predict the conductance for graphene junc- 
tions with leads made of other metals, and we present 
results for the conductance of junctions with Au leads. 

We perform transport calculations at zero source- 
drain bias using the nonequilibrium Green's function 
(NEGF) SMEAGOL code [l3], which is interfaced with 
the density-functional-theory (DFT) SIESTA package 
[l8| . We employ norm-conserving pseudopotentials [l^l 
in the local-density approximation (LDA) [l^l, and a 
real-space grid equivalent to an energy cutoff of 310 Ry. 
We explicitly construct fine-tuned basis sets for C and 
Al atoms following the prescription of Junquera et al. 
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The largest number of atoms (numerical orbitals) 
included in the present calculations is 464 (5600). We 
made modifications to the transport code to improve 
memory allocation and parallelization, in order to handle 
calculations of this scale. 

Graphene junctions without metal contacts. First we dis- 
cuss pristine graphene with no metallic leads attached. 
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FIG. 1: (a) Band structure and (b) density of states (DOS) 
for graphene near the Dirac point {Ed ~ 0). The two van 
Hove singularities are marked, (c) Conductance of a graphene 
ribbon 98.1 nm wide without metaUic leads. Inset: The resis- 
tance R diverges at the Dirac point. Dashed lines in (b) and 
(c) are results of tight-binding calculations, (d) Band struc- 
ture and (e) projected DOS on carbon atoms for graphene on 
Al. The highlighted lines in (d) are graphene bands. 



The band structure and density of states (DOS) in the 
vicinity of the Dirac point [Ed) are shown in Figs. [Ha) 
and (b), respectively. The DOS around Eu can be fitted 
to g{E) = D\E - Ed\ with D = 0.11/(eV2 unit cell). In 
all transport calculations described in this paper, an ef- 
fective ribbon width is determined hy W = nkW^, where 
Uk is the number of k-points used in the calculation and 
wq is the size of the unit cell in the transverse (x) direc- 
tion. We show in Fig. [Ijc) the conductance G as a func- 
tion of electron energy for a graphene ribbon oiW ^ 100 
nm. At this width, the quantized plateaus in the conduc- 
tance become unresolvable within the energy resolution 
employed (0.01 eV), and G resembles the linear behavior 
of the DOS around Ed. The boundary conditions and 
type of edge - zigzag or armchair ~ have an impact on 
the conductance only at the vicinity of the Dirac point, 
in a narrow energy range inversely proportional to the 
ribbon width. Away from Ed, the conductance remains 
unaffected by this choice with W being the same as the 
geometric width in the widest ribbons studied. Tight- 
binding results with a nearest-neighbor hopping param- 
eter 7 — —2.65 eV determined from the LDA bandstruc- 
ture are shown by dashed lines in Figs, [ijb) and (c) for 
comparison. The agreement with the DFT results is quite 
good within ±1 eV of the Dirac point. 
Junctions with Al-graphene contacts. The electronic 
structure of graphene in contact with Al (henceforth re- 
ferred to as the lead) shows a certain amount of hy- 
bridization arising from the interaction between graphene 
and the metal [Fig. Hl^d)]. The leads are modeled by 
six layers of Al in the (111) orientation with a graphene 
monolayer placed at a distance of 0.34 nm away, in the 
minimal-energy configuration [l6| . Given the relative po- 
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FIG. 2: (Color online) (a) Schematics of Al-graphene junc- 
tions. The semi-infinite leads, shown in the shaded regions, 
include both graphene and Al. The edges of the ribbon are of 
either armchair (1) or zigzag (1'). The unit cell sizes along 
the a;-direction are wq = 0.49 nm and w'q = 0.85 nm, respec- 
tively. L is the distance between Al atoms at opposite sides of 
the junction, (b) Energy location of the Dirac point with re- 
spect to the Fermi level across the junctions. Junctions 1, 2, 
and 3 have a fixed width W = 98.1 nm, and lengths L = 3.40 
nm, 6.80 nm, and 13.60 nm, respectively. The DFT results 
(dots) are shown along with fitted curves (see text). 



sition of the Dirac point Eo with respect to the Fermi 
level Ep at the lead, A — Ed — Ep ~ —0.6 eV, graphene 
becomes n-doped in the contact region. The DOS pro- 
jected on C atoms at the leads gc is nonzero at the Dirac 
point with a finite value of gc{ED) ^ 2 x 10^^ elec- 
trons/eV per graphene unit cell [Fig. [TJe)]. The fluctua- 
tions of gc above Ep are due to hybridization. 

The junctions studied are schematically shown in 
Fig. [21 The leads are displayed within the shaded ar- 
eas. The Hamiltonian, overlap, and density matrices of 
the leads are obtained self-consistently from DFT bulk- 
calculations and used to set up the electron density, 
chemical potential, and self-energies for the transport cal- 
culations. Since the leads include graphene, the contact 
area between metal and graphene is essentially infinite. 
This enhances the transparency [9|] , eliminating the large 
contact resistance that would arise if charge carriers were 
to tunnel from Al to graphene. An additional number of 
Al atoms between the leads and the freestanding section 
are added so that the self-consistent electrostatic poten- 
tial across the junction develops smoothly. We have in- 
vestigated various configurations to understand the ef- 
fects of the length (L), width (W), and orientation of 
the graphene ribbon on electron transport. Junctions 1, 
1', 2 and 3 have lengths L = 3.40 nm, 3.44 nm, 6.80 
nm, and 13.60 nm, respectively. The effect of varying 
L is studied through junctions 1, 2, and 3 with a fixed 
width oi W = 98.1 nm [see Fig. [2I^b) for geometries]; the 
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FIG. 3: (Color online) (a) Conductance of junctions with a 
width ofW = 98.1 nm and different lengths L. The curves are 
offset by 20 x 2e^ /h for clarity. Two prominent conductance 
minima at En (—0.6 eV) and Ef (set to zero) develop as L 
increases. The symbols represent first-principles results, and 
the solid lines are from model calculations (see text). Inset: 
minimum conductance versus L (the dashed line is a guide 
to the eye), (b) Conductance of junctions with a length of 
L = 3.4 nm and different widths W , highlighting the opening 
of a conduction gap as W decreases and a sharp peak pinned 
at Ed = -0.6 eV. 



conductance G as a function of W is studied for junction 
1 (L = 3.40 nm); and the effect of ribbon orientation 
(armchair versus zigzag edges) is studied for junctions 1 
and 1' by choosing a similar width for them: W = 98.1 
nm (1) and W — 98.6 nm (1'). (The potential profiles 
in Fig. [2Ib) will be discussed later.) 

The conductance for junctions 1-3 (with a fixed width 
of = 98.1 nm) is shown in Fig. [2Ja). One notices the 
lack of symmetry with respect to both the Fermi level 
(energy zero) and Ejj in the lead (—0.6 eV) in all cases 
regardless of the particular values of L, in contrast to 
the symmetric curve for pristine graphene in Fig. [TJc) 
and those displayed in Ref. 0]. Similar asymmetries 
in the conductance have been seen experimentally, e.g. 
Refs. [l^ and [2^. In addition to the conductance min- 
imum at Ed, the emergence of a second one at Ep is 
apparent, and it becomes more prominent as the length 
increases. This is because, as i ^> oo, the electron DOS 
approaches g{E) at the freestanding part of the junc- 
tion and the conductance becomes that of two resistors 
in series GiE) oc giE)gciE)/[giE) + gdE)], where 
gc{E) cx £>|£' — A| is the projected DOS of the clamped 
region. Based on the same underlying physics, a double 
peak in resistance was experimentally demonstrated for 
graphene devices comprising two regions with noticeably 
different doping levels [23]. The two conductance min- 



ima were also observed in junctions with gate-tunable 
barriers 2j| . The inset in Fig. [Sja) shows the minimal 
conductance near Ed for the values of L studied; it is of 
the order of 2e^//i. A small peak at Eo (—0.6 eV) is vis- 
ible for the two shorter junctions 1 and 2, and becomes 
unnoticeable for the longer junction 3. This small peak is 
due to states that exist only in the contact areas through 
hybridization with Al atoms and decreases rapidly with 
increasing junction length. 

In Fig. 131b) the conductance is plotted for the shortest 
junction [L = 3.4 nm) with four values of width W: 9.8 
nm, 24.5 nm, 49.1 nm, and 98.1 nm. The logarithmic 
scale of the conductance helps examine the linear depen- 
dence of the conductance on W for energies sufficiently 
away from En. As mentioned above, the peak at —0.6 
eV only appears in very short junctions. This is yet to 
be observed in experiment, where longer ribbons are cus- 
tomarily employed 
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Set aside this peak, the open- 
ing of a conduction gap with decreasing W is evident, 
as well as the significant and nonlinear drop in conduc- 
tance for the narrowest junction with W = 9.8 nm, in 
full accordance with experiment 25[. To investigate the 



effect of the edge orientation, we study the conductance 
of armchair (1) and zigzag (1') junctions of comparable 
dimensions [Fig. |4lja)]. Only small deviations are found 
resulting from the anisotropy of the band structure. The 
effect of orientation is quite small in the vicinity of En . 
Tight-binding model. We next present a tight-binding 
model that can be used to study junctions bigger in size 
and with arbitrary nonbonding metal contacts. To that 
end, we estimate the effective potential profile V{z) at the 
atomic scale for 7r-electrons along the junction, directly 
related to the spatial dependence of doping created by 
the metallic leads. This profile can be obtained semiclas- 
sically from DFT as the value of the energy at the Dirac 
point with respect to the Fermi level, as a function of z: 
V{z) — Eo^z) — Ep. We extract this quantity from the 
position-dependent energy shift of the van Hove singular- 
ities closest to the Dirac point [denoted by vHi and vH2 
in Fig.[IJb)] along the junction. The results for junctions 
1-3 are shown in Fig. [2Jb). Significant doping occurs at 
distances up to only a few nanometers from the edge of 
metallic leads. For the two longer junctions, 2 {L =6.80 
nm) and 3 {L =13.60 nm), Ep) is very close to Ep in the 
middle of the freestanding section {z ^ 0). We fitted our 
DFT results for V{z) by the following expression [solid 
lines in Fig. [2Ib)]: 



V{z) 



A cosh(z/A)/cosh(Lcff/2A) 
A 



for |z| < icff/2, 
otherwise. 

(1) 

The parameters used are A = —0.6 eV, A = 1.05 nm 
and Leff = L -\- 5ao, where ao = 0.142 nm is the C-C 
interatomic distance. We emphasize that on the sub- 
nanometer scale, V{z) starts to deviate from its asymp- 
totic value inside the region of intermediate Al atoms. 
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FIG. 4: (Color online) (a) Conductance for junctions 1 (arm- 
chair edges, L — 3.40 nm, W = 98.1 nm, light line) and 1' 
(zigzag edges, L — 3.44 nm, W — 98.6 nm, dark dots), (b) 
Conductance for junctions consisting of graphene with Au 
leads for two values of L from model calculations (see text). 
The junction width is = 98.1 nm. 



before reaching the atomic boundary of the contact in- 
dicated by the vertical dashed line in Fig. [2][b) . In our 
calculations, sufficiently many Al atoms are included be- 
tween the bulk (shaded) leads and the freestanding sec- 
tion in order to map the transitional behavior of V{z). 
This profile provides a refinement of models where sharp 
steps at the metal edge are assumed [1, and will be- 
come relevant as experimental junctions shrink in size. 

We have performed conductance calculations within a 
tight-binding approach using the potential profile V{z) 
and a smearing 6 = 8 meV introduced through the leads 
self-energies [8| (to account for the nonzero conductance 
at Ed and Ep and the small peak at Ed in Fig. 
We obtain excellent agreement with our full-scale calcu- 
lations [as indicated by the solid lines in Fig. [3l[a)], with 
the exception of sharp dips in the conductance [indicated 
by arrows in Fig. O^a)] which are not reproduced. These 
dips occur at the anticrossings energies due to hybridiza- 
tion with aluminum [cf. band structure in Fig.[TJd)]. We 
have also used the model to calculate the conductance for 
a longer junction {L = 40 nm), where the two conduc- 
tance minimum are fully developed. In the limit A — )■ oo 
only one conductance minimum appears and as a conse- 
quence the electron- hole symmetry is preserved [Hj. 

Junctions with other nonbonding metal contacts. We use 
the model presented above to make a prediction for the 
conductance with Au as the metal lead. In this case 
A = -1-0.12 eV and graphene is p-doped [l^. Using the 
same values of A, LeS in Eq. (1) and the smearing pa- 
rameter 5 as before, the conductance from tight-binding 
calculations is shown in Fig.m^b) for two junctions with a 
fixed width oiW = 98.1 nm. The solid line shows results 
for L = 13.6 nm and the dashed line for L = 40 nm. For 
small values of A, the two conductance minimum may 
be smeared out if the amount of charge fiuctuations is 
significant 



12l . Il5l . |26|. Nevertheless, the electron- hole 
asymmetry should remain noticeable in the conductance 
curve. Our results give microscopic justification to prior 
theoretical studies where graphene junctions are modeled 



as an isolated graphene sheet with a space-dependent po- 
tential 0. However, we have demonstrated that to ob- 
tain accurate results, one needs to take into account a 
small energy broadening [Si and allow for a finite poten- 
tial at the contacts. 

Conclusion. We have performed transport calcula- 
tions for graphene junctions attached to Al leads using 
first-principles nonequilibrium Green's function methods. 
The conductance features vary with the length and width 
of the junction, but are less sensitive to the ribbon orien- 
tation. We show that nonbonding metallic leads induce a 
noticeable electron-hole asymmetry in the conductance. 
The opening of a conduction gap with decreasing W is 
also captured. Two minima in the conductance emerge 
for large enough junctions. In addition, our calculations 
yield reliable information on the doping variation along 
the junction for metallic leads interacting weakly with 
graphene, and we find accurate potential profiles at the 
vicinity of the metal boundary. We demonstrate that the 
dominant features of our first-principles results can be 
reproduced by an analytically tractable effective model. 
The parameters of the effective model are derived from 
first-principles. As an application we use the effective 
model to predict the conductance in junctions with Au 
as the supporting metal at the leads. 
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